Loop Quantum Dynamics of the Schwarzschild Interior 
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Wc examine the Schwarzschild interior of a black hole, incorporating quantum gravitational mod- 
ifications due to loop quantum gravity. We consider an improved loop quantization using techniques 
that have proven successful in loop quantum cosmology. The central Schwarzschild singularity is 
resolved and the implications for the fate of an in-falling test particle in the interior region is dis- 
cussed. The singularity is replaced by a Nariai type Universe. We discuss the resulting conformal 
diagram, providing a clear geometrical interpretation of the quantum effects. 
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O ' I- INTRODUCTION 



A long held expectation is that the singularities predicted from general relativity signal a breakdown of the classical 
theory requiring a more proper accounting of the quantum effects of gravity. The situation is analogous to Hydrogen 
q-i atom in ordinary mechanics, where classically the electron is expected to in-spiral towards the proton leading to a 
singularity. This behavior is cured when quantum mechanics is taken into account leading to a stable non-singular 
Hydrogen atom. In general relativity, two of the most relevant forms of singularities are the big-bang cosmological 
singularity, as well as black-hole singularities. A central question to be answered by any quantum theory of gravity 
is whether these singularities are regulated by incorporating quantum effects. Furthermore, if a quantum theory 
of gravity were to resolve the classical singularities, it would be of paramount interest as to what replaces them 
■ dynamically. 

\ A leading candidate theory of quantum gravity is known as loop quantum gravity (for reviews see [1-3]). The 
application of loop quantum gravity techniques to cosmological models is known as loop quantum cosmology [4] 
which has led to a resolution of the big-bang singularity, replacing it with a big-bounce for homogeneous and isotropic 
, models [5-7] . These results provide tantalizing first hints that loop quantum gravity may indeed provide the singularity 
resolution hoped for in a quantum theory of gravity. Thus the next step is to consider loop quantum gravity for the 
, black hole scenario. 

The simplest first step in considering loop quantum black holes consists of examining the interior of a Schwarzschild 
black hole. There, the temporal and radial coordinates flip roles, and the interior becomes a Kantowski-Sachs cosmo- 
logical model whereby the metric components are homogeneous and only time dependent. The interior therefore can 
be quantized in a similar fashion as loop quantum cosmology leading to the possibility that the singularity is resolved 
as in the cosmological case. The loop quantization of the Schwarzschild interior has been initially developed in [8] 
and [9]. There, the quantization indicates that the quantum Einstein equations are non-singular in a similar way loop 
quantum cosmology was originally shown to be non-singular [10]. However, the question of what replaces the black 
hole singularity is not answered. 

Therefore, we attempt to provide an answer as to what dynamics are predicted from the loop quantization of the 
Schwarzschild interior. Initial attempts in this direction are given in [11] which is based on the quantum Einstein 
equations of the original quantization of [8]. However, this quantization is the direct analog of the original formulations 
of loop quantum cosmology which have been shown to not have a good semi-classical limit [5]. In particular, a 
crucial parameter in the construction of the Hamiltonian constraint operator is taken to be a constant in the early 
constructions and it is this assumption that can spoil the semi-classical limit in loop quantum cosmology. An improved 
quantization of loop quantum cosmology whereby the parameter is assumed to explicitly depend on the scale factor 
has been applied in [5] and leads to the correct semi-classical behavior as well as more physically sensible quantum 
corrections to the standard cosmological dynamics. 
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In this paper, wc consider the improved quantization scheme of loop quantum cosmology applied to the loop 
quantization of the Schwarzschild interior. We will consider the resulting modifications to the classical interior by 
examining effective semi-classical equations that incorporate quantum effects. While we do not have a complete 
handling of the physical sector of the quantum theory owing to its complexity, the effective semi-classical description 
is intended to be an approximate shortcut to the physical predictions of the quantum theory. In the framework of the 
effective description, we will show that the singularity is resolved and discuss the implications for an in-falling test 
particle as well as discuss the resulting conformal diagram for the non-singular interior region. 



II. CLASSICAL FRAMEWORK 



A. Classical Hamiltonian 



The interior of the Schwarzschild geometry can be described by a Kantowski-Sachs homogeneous model whereby 
the line element is given by 



ds 2 = -N(t) 2 dt 2 + g xx (t) dx 2 + gnn(t) dVL 2 



(1) 



with N(t) being the freely specifiable lapse function and dfl 2 representing the unit two-sphere metric given in polar 
coordinates 



dn 2 = de 2 + sin 2 e < 



(2) 



The topology of the spatial slices is R x S 2 hence the coordinates range from x £ R, 9 £ [0,7r], <f> £ [0, 2ir], In the 
standard Schwarzschild solution, the metric components are given by 



N(t) 2 



2m 



- 1 



gxx{t) 



2m 



- 1 



gnn{t) = t 2 



(3) 



for t in the range t £ [0, 2m], where m is the mass of the black hole. 

Loop quantum gravity is based on a connection-triad Hamiltonian formulation of general relativity Considering 
connections and triads that are symmetric under the Kantowski-Sachs symmetry group R x SO(3), the connection 
A\ and triad Ef after gauge fixing of the Gauss constraint are of the form [8] 



A\ Tl dx a 

E?nd a 



cr^dx + bT-2d0 — br\ sin 
p c T3 sin d x + pbT2 sin 9 



f T3 cos 



(4) 
(5) 



where the dynamical phase space variables b,c,pb,p c are all only functions of time, and Tj are SU(2) generators 
satisfying [n, Tj] — ^tjTk- The relationship between the triad variables pb,p c and the metric variables g xx , gan is given 
through the general relation 

which cxplictly implies 



E?E? = det(g) q ab 



g 

yxx — |~| 

\Pc\ 

gnn = \p c \ 



(6) 

(7) 
(8) 



We thus see that the triad component p c directly determines the physical radius of the two sphere which is proportional 
t° \/gnn- We will use this fact when interpreting the quantum dynamics. 

The classical Hamiltonian consists of a single constraint known as the Hamiltonian constraint. The explicit form 
of the Hamiltonian is given by 



H 



-N 



8ttG 7 2 



2bdVWc + (b 2 +l 2 )^ 



dx dtt . 



(9) 



Performing the spatial integrations we run into a problem as the integral over x diverges. Following [8] , if wc restrict 
the integration over a; to a finite interval Lq, the Hamiltonian becomes 



H 



-N 



2G 7 2 



2bc^p- c + (b 2 +1 2 )^ 



(10) 
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written in terms of rescaled variables 

b = b, c= L c, p b = L p bl p c = p c (11) 

The Poisson structure of these phase space variables is given by the only non- vanishing Poisson brackets 

{b, Pb } = G 1 (12) 
{c, Pc } = 2G 1 (13) 

Written in terms of these variables, the metric (1) is given by 



ds 2 = -N 2 dt 2 + ^-^ T + \p c \dn 2 (14) 
\Pc\ J^o 



The equations of motion are given as in an ordinary Hamiltonian dynamical system based on the Hamiltonian (10). 
The dynamical behavior of the phase space variables is governed by Hamilton's equations which give for instance 

p b = { Pb ,H} = -G 1 — (15) 

and so forth. Additionally, the lapse -/V appears in the Hamiltonian as a Lagrange multiplier and this implies that 
the Hamiltonian must vanish leading to an additional constraint on the dynamics 

H = 0. (16) 

We will next consider two classes of solutions that will be relevant to the discussion. The first consists of the usual 
Schwarzschild interior we are most interested in. The second corresponds to the Nariai solution which will be of 
relevance when we discuss the modified quantum dynamics. 

B. Schwarzschild solution 

To arrive at the Schwarzschild interior solution let us first choose for simplicity the lapse (for b ^ 0) to be given by 

N = 2y^i (17) 



whence the Hamiltonian becomes 



The equations of motion are given by 



H 



2G 7 



2cp c + {b + 1 2 /b) Pb . (18) 



b = -I(6 + 7 2 /6) (19) 

c = -2c (20) 

P b = f(l-7 2 A 2 ) (21) 

Pc = 2pc (22) 

H = 0^2c=-^(b + ~f 2 /b). (23) 



The solution to these equations of motion is given by 



Pb(T) = p { b 0) e T V2me~ T - 1 (24) 
Pc(T) = e 2T (25) 



b(T) = ±7V2me- T - 1 (26) 
c(T) = - 7 mp[ 0) e- 2T . (27) 
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with ' a constant of motion. Transforming to a new time coordinate T = \ogt, we have NdT — 7-^p ^ = 
1/ yj2m/t — 1 dt and the metric (14) is given by 

1 / (0) \ 2 „ 

ds 2 = - 1 ^— [ de + - \)dx 2 +<W (28) 

If we then identify = £0, wc arrive at the Schwarzschild interior solution (3). 

In these coordinates, the singularity is located at t = 0, while the horizon is at t = 2m. We can identify the 
singularity and horizon in a coordinate invariant fashion by noting that at the singularity both triad components p c 
and pb vanish. The horizon then occurs where pb — and p c = Am 2 . 

Before moving to the second class of solutions let us discuss certain gauge freedoms of the Schwarzschild solution. 
Classically there are two gauge freedoms corresponding to coordinate changes in either t or x variables. A change in 
the time coordinate t leads to a change in the lapse leaving all physics unchanged. On the other hand if we consider a 
coordinate rcscaling of x, L n is correspondingly rcscaled such that the ratio dx/L is left invariant. This implies that 
a rescaling of x does not rescale the value of pb- Note that this point was mistaken in [8] where it was mentioned that 
a rescaling of p h ^ can be accommodated by a rescaling in x. 

In addition to the classical coordinate gauge freedoms, there is an additional gauge in the loop setup pertaining 
to the choice in restriction of the spatial integration leading to the Lq parameter. From the definitions of pb and c 
in (11), it is easy to see that if we choose a different size of Lq then this will rescale the value of pb and c and hence 
in terms of the Schwarzschild solution, p^ will rescale accordingly. 

The loop quantization is based on the phase-space variables pb,p c ,b,c all of which are invariant under a coordinate 
rescaling of x. The quantization is thus manifestly invariant under the first gauge transformation. However the 
quantum effects in principle can depend on the second gauge freedom which does not appear classically. This will be 
manifested in the quantum dynamics by a dependence of the physical results on the parameter p^ . Any dependence 
of the quantum physics on this gauge freedom implies the need for additional input into the theory (namely the 

specification of the value of pP' ) which must be done by hand. We will show that the improved quantization we focus 
on in section IIIB will be independent of and hence will not require additional input of its value. 



C. Nariai solution 



An additional classical solution that will be of relevance to the discussion of the effective dynamics is the Nariai 
solution [12, 13]. It differs from the Schwarzschild solution in that the two-sphere radius determined from p c is 
assumed to be constant and the model is sourced by a cosmological constant with the possible inclusion of charge. 
The assumption that p c is a constant requires that 6 = and hence the lapse (17) used in the Schwarzschild case 
becomes degenerate. We must therefore choose another form of the lapse to generate the Nariai solution. 

Choosing the lapse to be equal to one and adding in a cosmological constant A leads to the Hamiltonian 



H = 



-1 



2G 7 2 



Pb 



2ky^+(& 2 +7')^= 



/Pc 



A 
2G 



PbVPc- 



(29) 



We are interested in a solution where the radius of the two-sphere is constant implying p c = const, such that p c = 0. 
The resulting equations of motion imply 6 = 0. Putting 6 = into the equation for 6 yields the relation between 
the constant radius of the sphere and the cosmological constant p c = 1/A. The remaining two equations of motion 
simplify to 

(30) 
(31) 



c = 7AV Apb 

1 1 

Pb = — 7=c 
7 VA 

which can be written as a single second order differential equation for pb given by 

p b = Ap b . 

The general solution is given by 



(32) 



Pb = c\ cosh(v / At) + C2 sinh(v / At) 



(33) 
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where ci and C2 are two constants of integration. To obtain the global cosmological metric of the Nariai Universe we 
set ci = 1/ (A) and c 2 = to get 

p b = icosh(\/Ai). (34) 
Hence, we obtained the following solution of the classical field equations 

ds 2 = -dt 2 + — cosh 2 (VAt)^ + T dn 2 (35) 
A Liq TV 

compare for example [14, 15], which corresponds to the uncharged Nariai Universe. 

Geometrically Nariai type spacetimes can be understood as four dimensional submanifolds of a flat six dimensional 
Lorentzian manifold being the product of two (three dimensional) spaces of constant curvature. If the respective 
length scales are denoted by a and b say, the generic Nariai type metric takes the form 

ds 2 = a 2 (-dT 2 + cosh 2 Td X 2 ) + b 2 dQ 2 . (36) 

For a 2 = b 2 = 1/A, this corresponds to the uncharged Nariai universe, given by (35) upon a rescaling of the time 
t = t/a. If a 7^ b, the spacetime corresponds to the charged Nariai universe. It should be noted that the knowledge 
of the metric alone, does not suffice to distinguish between the two types of possible charges, electric or magnetic. 



D. Radial geodesies 

In order to interpret the quantum dynamics, we will also consider the trajectory of a radially in-falling test particle. 
The geodesies can be derived from the line element 

ds 2 = -N(t) 2 dt 2 + y~~7777 + \Pc(t)\dn 2 . (37) 

Since we are mainly interested in radial geodesies we henceforth assume dfl 2 = constant. The geodesic equations are 
most easily read off from the Lagrangian 

dt\ 2 p 2 ( dx^ 



where C = for a masslcss particle and C = 1/2 for a massive particle, and r refers to the proper time for a massive 
particle. Since the Lagrangian (38) is independent of the variable x, its conjugate momentum tt x is conserved 

n x = — = p2 = const = -8 (39) 

ax dr PcLq 

where £ represents the total energy including the gravitational potential energy for a static observer at infinity. Since 
dx/dr = £p c Lq/pI the Lagrangian can now be written as 

. 2 C = -N 2 (^-) 2 + ^£ 2 . (40) 



dr 



PI 



Hence, we obtain 

In the classical Schwarzschild case, we have p c L\j (N 2 p\) = 1 and hence the geodesic equation can be understood 
as a point particle with energy £ moving in an effective potential given by 2C/N 2 . Moreover, the geodesic equation 
is simple to interpret because of the simple relation between the coordinate parameter t and the two-sphere radius 
determined fromp c since we havep c = t 2 . However, when we come to the quantum dynamics, the relationship between 
the coordinate parameter t and p c will become complicated and we will no longer be interested in the behavior of 
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the coordinate parameter t as a function of proper time. We will therefore be most interested in p c (r) which can be 
derived from 

dp c = dpc_dt_ , 
dr ~ dt dr [ ' 

where dp c /dt is to be determined from the equations of motion derived from the Hamiltonian, and dt/dr is given 
by (41). 

Classically, for any given energy £, the in-falling particle reaches the singularity at p c = in finite proper time. 
We will be interested in whether the quantum dynamics implies that the singularity is never reached by the in-falling 
particle. We will see that the loop quantum dynamics imply a minimum value of the two-sphere radius p c and hence 
the solutions are non-singular with an in-falling particle never reaching the classical singularity. 



III. EFFECTIVE LOOP QUANTUM DYNAMICS 



In attempting to ascertain the quantum corrections to the classical behavior implied by the loop quantization, we 
will work in an effective semi-classical description that leads to modifications of the classical Hamiltonian constraint. 
The effective description is motivated from the construction of the quantum Hamiltonian operator, although a more 
rigorous understanding of the quantum dynamics would require more sophisticated machinery. Our aim here is to 
describe the expected dominant corrections arising from the quantum theory and as such, there is room for additional 
corrections that could modify some of the results presented here. Evidence that the effective theory we describe 
provides an accurate description of the quantum dynamics is given by homogeneous and isotropic loop quantum 
cosmology with a massless scalar field with or without a cosmological constant [5-7] . 

The main quantum correction arises from the fact that in the loop quantization no direct operators exist corre- 
sponding to the connection components b and c. Instead, the quantization proceeds as in loop quantum gravity by 
considering holonomics of the connection. The holonomics roughly consist of exponential terms like cxp(ib5b) with 
5b corresponding to the edge length of the holonomy. The end result in the effective Hamiltonian is to replace the 
classical b terms in the Hamiltonian (10), with sin(b5j,)/Sb terms for instance. Thus the 8b, 5 C factors play a crucial 
role in determining the effects of the quantum corrections. The classical limit is recovered by taking the Sb,S c — > 
limits. The Sb,S c parameters, arising from the holonomy edge lengths, are a measure of the quantum discreteness 
arising from the loop quantization. 

In the original loop quantization of the Schwarzschild interior [8], the Sb,5 c terms were assumed to be a constant 
value labeled by 6 analogous to the original quantization of loop quantum cosmology where similar terms appear in 
the Hamiltonian operator. In [5], it was shown that because of the assumption that 5 is a constant (in loop quantum 
cosmology the parameter is labeled ^0)7 loop quantum cosmology does not have the correct semi-classical limit. The 
solution to this issue presented in that article is to relax the assumption that the S factors are constant and instead 
consider them to be functions of the triad variables. The precise prescription given leads to the non-constant factor 
labeled \x which varies inversely with the scale factor of the universe and with this behavior the correct semi-classical 
limit is achieved [5]. 

Thus our goal will be to consider the effective dynamics with a quantization scheme arising from the improved 
prescription for loop quantum cosmology given in [5]. In particular, the Sb,S c factors will depend explicitly on the 
triad components pb and p c . We will start however, with the case of constant 5 since the effective dynamics are exactly 
solvable and allows to contrast the modifications from those arising from the new non-constant S quantization. The 
results for the constant 5 case have been previously derived in [11]. We merely repeat the results that are relevant to 
our discussion. 



A. Constant S Quantum Hamiltonian 



With the considerations of the previous discussion the quantum modifications can be understood as leading to 
an effective Hamiltonian with the holonomy parameter 6 playing a crucial role in determining the magnitude of the 
corrections. The effective Hamiltonian constraint reads [11] 



H 5 



N 



2G 7 2 



sin b5 sin cS 



sin 2 M 



Pb 

\/P~c 



(43) 



It is easy to see that in the limit 6 — ► 0, the classical Hamiltonian (10) is recovered. In analogy to the classical case, 
we choose ./V = j^/p^ : 5/ (sin b5) and the equations for the pairs (b,pb) and (c,p c ) decouple to two independent pairs of 
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differential equations 



c=2 7 G^ = -2 — (44) 
OH 

p c = -2 7 G-^ = 2p c cos c5 (45) 
• dH s f[\ 

& = 7 ^ = -2(^ + ^mJ (46) 

dH s 1 , f /, 7 2 <5 2 \ 

Pb = -1st = o cosW 1 — nrr Pt' ( 47 ) 



96 2 V sin 2 

The equation of motion for c can be integrated to give 

c(T) = i arctan(=F7TOpi 0) ^ e_2T ) ( 48 ) 

where the constant of integration has been fixed so that the limit of the effective c(T) as 6 — > coincides with the 
classical expression (27). Next, the equations for p c can be solved 

Pc (r),e 2T + 7 V(pf»)^ e - 2T (49) 
Also the equation for b can be integrated and yields 

( . / T* — 1 no' 9m r^/i\ 1 

(50) 



6(T) = - arccos 
o 



v/TTWtanhj ^( ^'f'" - log ^) } 



Finally, Pb(T) can be obtained from the vanishing of the Hamiltonian (43) which gives 

sin b5 sin c5 2p c 
Pb = x x ■ 277775 I — a- (- 51 ) 

/From the form of 6(T), it is evident that the evolution ends at a minimum value T m j n where the absolute value of 
the argument of the arccos function equals one. Moreover, from the form of p c (T) it is evident that the group orbits 
also have a minimal value at T = T m j n = log( r ymp^6/2)/2, with p c reaching its maximum value at the classical 
horizon and at T m - m . This minimal two-sphere radius is given from p c (T m - nl ) = jmpf 5, which depends on the mass 

of the black hole and also on p£ . Figure 1 shows the plots of pi, = Pb{p c ) for various values of m and p^ . The 
singularity is thus avoided through a bounce in the two-sphere radius and the quantum dynamics matches two black 
holes together though in a generically non-symmetric value. A symmetric solution can be constructed by fine tuning 
the value of p£ ■ 

As we have seen in the figures, the mass parameter m controls the mass of the first black hole whereas p^ controls 

the mass of the second black hole. This suggests that in the effective quantum theory, p^ would not be regarded as a 
meaningless constant of integration or a certain choice of gauge. For the constant 5 quantization of the Schwarzschild 
interior spacetime this parameter acquires an important physical meaning, and two different choices correspond to 
different effective quantum theories. Any physical results derived from the constant S quantization must include a 
prescription for specifying the value of p£ which has to be input by hand. 



B. Improved Quantum Hamiltonian 

We now wish to relax the assumption that the 5 parameters are constant. In accordance with the results from loop 
quantum cosmology [5], we will assume that 8b, S c depend on the triad variables p c and pb- The method to constrain 
the exact dependence relies on the fact that the holonomies used to construct the quantum Hamiltonian operator 
form closed loops. The behavior of 5b and S c can then be constrained if we demand that the physical area of the 
closed loop be equal to the minimum area gap predicted by loop quantum gravity 



A min = A = 2 V / 37T7 $ 



(52) 
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FIG. 1: These figures show the effective plots of pj, = Pb(Pc) for the case with 8 = const. The dashed lines represent the 
classical solution while the solid lines correspond to the solution of the effective quantum theory with S = 0.1. Clearly, the 
initial singularity is removed and the spacetime contains a second black hole. The three plots on the left show the solution for 
varying mass parameter, m = 1; 1.5; 2 respectively. On the right p^ ' is varied in the three plots, = 1; 2/3; 1/3 respectively. 
One clearly notes (left) that the mass parameter m controls the mass of the first (classical) black only, while (right) the 
parameter p^ controls the mass of the second black hole. 



This procedure can be executed in different possible ways in the Kantowski-Sachs model considered here, whereas 
homogeneous isotropic models lead to a unique choice. We discuss in this section what we feel to be the most 
geometrically natural scheme. We will show in the section, that this scheme has the advantage that the resulting 
dynamics is independent of the p£ parameter. We will discuss an alternative scheme in appendix A. 

The procedure considered in this section amounts to evaluating the classical area of the holonomy loops and 
constraining the area based on this. For instance if we consider the holonomy loop in the x, 9 plane, the classical 
physical area is given by 



SbS c Pb 



(53) 



which is evident from the form of the metric (14). For the holonomies on the two-sphere, the loop does not close as 
discussed in [8] , but we can assign an effective area given by 



If we constrain these areas to be equal to A„ 



Adc/> = $bPc ■ 

we have the following behavior 



5 b 

So 

With this, the effective Hamiltonian constraint reads 

N 



A 



Pb 



H s = 



2G 7 2 



sin b6(, sin 5 c c 

2 — J j — VPc 



'Sin 



7 



Pb 

\/P~c 



(54) 

(55) 
(56) 

(57) 



Once again, in the limit A — > the effective Hamiltonian (57) reduces to the classical Hamiltonian (10). 

We can now derive the equations of motion as in the previous case. Choosing the lapse function to be N = 
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7 v /p2>fc/(sinM;>) ! the resulting equations of motion are given by 

sin c8 c 1 sin Sbb 



- ccoscSc 

o c I o c 



+5^«»w^-^-) + Tf-V (58) 

2 o c V sin 65b / 2 o c sin &Ob 

p c = 2p c cos c<5 c (59) 

1 /sin 65ft 7 2 5 b \ sinc<5 c 5 C 

-- — = — + 5 — +^-ccosc<5 c (60) 

2 V db sin 60{, / Ob Of, 

„2i 



1 ^ / 7 2 5 



6 

Pfa = ^cos6<5 6 (l , ' 2 " W (61) 

2 v sin 6o&/ 

Before considering numerical solutions to these equations, we can get an idea on where the quantum modifications 
should be appreciable. It is evident from the effective Hamiltonian, the the modifications arise when either c8 c or 
bSb are of the order of one. Near the classical singularity we have pb — > 0,p c — > 0, 6 — > ±oo and c — > — oo in such a 
way that 6<5b and cS c blow up, hence we expect the quantum corrections to become important near the singularity. If 
we also consider the horizon, we have pb — > 0,p c — ► (2m) 2 , 6 — * ±0 and c — ► -7mp[°'/(2m) 2 such that 6<5b is small, 
but co" c diverges. Thus we can expect quantum corrections also as the classical horizon is approached. Therefore, 
the picture that arises is a region in the interior where classical behavior is recovered with quantum effects near the 
classical singularity and also near the classical horizon. This behavior is borne out in the numerical simulations which 
we now discuss. We focus in this article on the effects near the singularity and will comment later on the horizon 
effects. Accordingly, in the numerical simulations we will specify initial conditions starting in the interior region where 
the quantum effects are negligible and evolve towards the singularity. 

To solve the equations of motion numerically, we start at some initial point near the horizon with initial conditions 
p c {To) , Pb{To) , 6(To) that match the classical solutions (25,24,26) for some initial To and the constraint H$ = is then 
solved to get c{Tq). The coupled differential equations of motion are then solved for the approach towards the classical 
singularity (T = -co). The results of the simulations for pb as a function of p c appear in figure 2 for various initial 

conditions corresponding to different values of the classical black hole mass m as well as different values of p b °^ . The 
results all indicate that the classical singularity is resolved with the two-sphere radius p c being bounded from below. 
The second figure of 2 indicates that choosing different values of only rescale Pb{T) as in the classical case, and 
thus the effective dynamics predicted is insensitive to this classical gauge choice as opposed to the constant 5 case. 

Unlike the constant 8 case of section III A, the solutions do not match up to a separate black hole solution. Instead, 
the solutions asymptote to a Nariai type solution for T — > — oo. The solution is characterized by 



6 = 6, Pb = Pbe 



-aT 



,-oT 



Pc = Pc, (62) 



where the barred quantities and a arc constants. Figure 3 shows this behavior with plots of p c , b and the ratio of 
c/pb as functions of time. It is clear that p c undergoes damped oscillations settling in to a fixed finite value as T goes 
to minus infinity. Similarly the ratio of c/pb settles into a finite value given by c/pb as expected from a Nariai type 
solution. 

^From the equations of motion we can determine the asymptotic values of the constants p c , 6, a, and c/pb- Assuming 
Pc = Pc is a constant, Eq. (59) implies 



cos c8 r = cos 



(Va^c) = 0. (63) 
V Pb ' 



cos (VA#c) =0, 1 = (64) 
V pb J Pb 2^/Ap c 

which also yield sin o~ c c = — 1. Assuming furthermore that 6 = 6 is constant, the resulting equation of motion for 6, 
namely Eq. (60) yields 

2sinfVAA) -Binfvft-L) 2 = (65) 
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FIG. 2: These figures show the effective plots of pt = Pb(Pc) for the case where Sb,S c are not constants. The dashed lines 
represent the classical solution while the solid lines correspond to the solution of the effective quantum theory. The initial 
singularity is removed, however, not with a second black hole. The three plots on the left show the solution for varying mass 
parameter, m = 1; 1.5; 2 respectively. On the right pj; ' is varied in the three plots, p^ ' = 1; 2/3; 1/3 respectively. One verifies 
that this quantization scheme does not depend on the choice of p b °^ and thus the effective theory is independent of this gauge 
freedom. 




-25 -20 -15 -10 -5 -25 -20 -15 -10 -5 



T T 

FIG. 3: In the left figure b(T) (solid) and p c (T) (dashed) are plotted. Both tends towards constants as T — » — oo. The right 
figure show the ratio c(T)/pb(T) also tends to a constant asymptotically. This asymptotic behavior is indicative of a Nariai 
type metric of the form given in equation (36). 



The same equation is obtained by evaluating the vanishing of the Hamiltonian constraint (57). Next, Eq. (61) implies 
that Pb(T) = pbe~ aT and fixes the constant a to be 




Finally, the above equations Eqs. (64)-(66) can be inserted in the equation of motion for c leading to 




For given value of A this equation can be solved for the value of the ratio bj ' ^pp~ c . A find root algorithm is necessary 
since an analytical solution cannot be obtained. The ratio b/^/p^ then specifics the constant in the exponent a by 
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Eq. (66). Moreover, with 7, A and bj yfp~ c given, the relation (65) can be used to extract the value of p c which in turn 
yields the value of b. Finally, Eq. (64) fixes the ratio c/pb- However, there is a remaining scaling freedom in c and pb 
since only their ratio is specified. Note that this implies that rescalings in pb have to be compensated by simultaneous 
rescalings in c. 

^From the previous equations it is clear that the asymptotic value of p c depends only on the constants A and 7 as 



Pc 



7 2 ff(A) 



(68) 



where g(A) is some function determined from the solutions of equations (67) and (65). For the natural choice 
A = 2v / 37T7Zp 1 the following asymptotic values are obtained 



b 

c/pb 



0.156, 

-2.290 mf,!, 



Pc 

N 



0.182/1,!, 
0.689, 



aw 0.670 



(69) 



where A" is the asymptotic value of the lapse which also behaves as a constant. These values agree with the asymptotic 
region of the numerical solution that we studied, see in particular figure 3. 





FIG. 4: Left: p c as a function of proper time of radially in-falling geodesic. The black dashed line represents the classical 
behavior indicating that the particle reaches the singularity in finite proper time. The solid line is the geodesic inside the 
quantum black hole. The particle approaches the singularity but undergoes damped oscillations at a fixed two-sphere radius 
eventually settling in at a fixed value p c . For the natural values of the parameters 7, A, the two-sphere radius is on the order of 
the Planck length. Right: The logarithm of the Kretschmann invariant as a function of the time coordinate T is plotted. The 
dashed line is the classical value which behaves as exp(— 6T) and diverges at the singularity T — —00. The solid line represents 
the quantum value which asymptotes to a fixed value which again indicates a resolution of the singularity. The numerical 



asymptotic values obtained for p c and the Kretschmann invariant K are p c 
independent of the black hole mass. 



0.182 Z?,, and K « 124.35 mi, which are both 



The fact that p c asymptotes to a fixed value on the order of the Planck length squared implies that the radius of 
the two-sphere asymptotes to a value on the order of the Planck length. Furthermore this value is independent of 
the mass of the black hole. The interpretation of the effective dynamics becomes clear if we consider the radially 
in- falling test particle using equation (42). A plot of the numerical integration of p c (r) is shown in figure 4. An 
in-falling particle classically would reach the singularity at p c = in finite proper time as evidenced by the solid black 
line, however as a result of the quantum effects, the particle oscillates about the fixed value of p c . The particle thus 
settles in at a fixed Planckian radius. The geodesic reaches this asymptotic value in infinite proper time and hence we 
can say that the singularity is truly absent in the quantum modified spacetime. The fate of an in-falling test particle 
is also of interest in the context of the information loss paradox. The results are suggestive that the black hole can 
store information infinitely long, and strictly speaking no information can ever be lost. However, once the thermal 
black hole radiation is taken into account, things are still more involved. In this respect, however, we would like to 
point out that the improved quantization scheme also leads to quantum modification of the horizon. In the context 
of Euclidean quantum gravity, it has been argued in [20] that indeed horizon modification might resolve the paradox. 

We have also considered the behavior of the Kretschmann invariant, the square of the Riemann tensor K = 
R a bcdR abcd - For the metric (14), written out explicitly, this yield a rather complicated expression. However, plotting 
this quantity for our numerical solution also shows that the geometry of the spacetime changes significantly, from 
the classical Schwarzschild like behavior where K = 48m 2 exp(— 6T) (which diverges at T — —00) towards a space of 
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constant curvature, as expected for the Nariai type metrics whose form is given by 

^ = 4 (^ + m)- (70) 
Plugging in the theoretical values in (69) gives an expected value 

Ktt 124.36 jt4i (71) 

which is in good agreement with the numerical results. Furthermore this asymptotic value is independent of the 
black hole mass. We find that the quantum effects become appreciable when the Kretschmann curvature invariant 
approaches the Planck scale. Figure 4 shows a plot indicating this behavior. 

The quantum effects removing the singularity can best be visualized by Carter-Penrose diagrams, see e.g. [16]. The 
classical diagrams for both the Schwarzschild and Nariai spacetimes are given in figure 5. The quantum effects remove 
the classical Schwarzschild singularity shown as the wavy line at the top of region II in the Carter-Penrose diagram 
and glue the diagram to the shaded region of the Nariai diagram. The resulting diagram of the interior is shown in 
figure 6. 




FIG. 5: The left figure shows the Carter-Penrose diagram for the Schwarzschild spacetime. J? and J?~ denote future and past 
null infinity, respectively. The usual Schwarzschild metric in curvature coordinates only covers region I. The Kantowski-Sachs 
form of the line element covers region II. The dotted lines represent the r — constant lines in the respective regions. The right 
figure shows the Carter-Penrose diagram of the non-singular charged Nariai universe. The x = and x = 2-7T lines are identified 
and the dotted lines stand for the constant two-sphere radii. The solid and dashed diagonal lines represent the disconnected 
future and past event horizon, respectively. Every point of this diagram represents a sphere of the same area. For later purpose 
we have indicated a shaded region. 



T = -oo 




FIG. 6: In region II we have indicated a specific r — constant line. This lines separates the region with decreasing group orbits 
(lis) from the region where the 2-spheres have (effectively) constant volume (I In), which defines the Nariai spacetime. 



IV. DISCUSSION 

At the level of the effective description we have presented, we have been able to answer the questions posed at 
the beginning, namely what replaces the singularity in the Schwarzschild interior. In all of the quantization schemes 
we have discussed, singularity avoidance arises naturally thus providing further evidence that loop quantum gravity 
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successfully resolves the classical black hole singularity. We have seen two possible replacements for the singularity, 
one in which quantum effects connect two separate (non-identical) black holes resembling a wormhole-like solution. 
The other solution predicting repulsive gravitational effects leading to an equilibrium two-sphere radius to which an 
in-falling particle would asymptote. 

The results are suggestive but can not be considered as a proof that the LQG resolves the black hole singularity. 
The effective equations considered here are not rigorously derived from the quantum theory and thus can only be 
considered approximations of the quantum dynamics. Since understanding the quantum dynamics is difficult, we do 
not have a complete handle on when the approximation might break down. As was mentioned however, the effective 
equations have approximated the quantum dynamics in the models where both can be developed completely, namely 
the homogeneous and isotropic models of LQC sourced with a massless scalar field. 

Some initial attempts at constructing semi-classical states for the constant 5 quantization have be made in [21] and 
[22] and it is worthwhile to compare our results with the previous work. In the works cited, numerical instabilities of 
the difference equation are found which bring into question whether any effective picture is valid near the singularity. 
In [21] it is argued that no proper semi-classical states exist unless the Immirzi parameter 7 is less the the accepted 
value in the LQG literature raising a further question. In this article, our aim has been to probe the improved, 
non-constant S quantization and thus we do not see a direct contradiction. In particular, the bound on the Immirzi 
parameter depends explicitly on the constant value of 5 and so it is not clear if a bound exists in the improved 
quantization. In addition, both papers evolve Gaussian semi-classical states with the quantum difference equation 
which seem to imply that semi-classicality breaks down near the singularity. The numerical results are based on using 
the triad component p c as an internal clock, however, our results suggest that p c bounces and hence would not play 
the role of a good clock. Thus to reconcile our results with those of [21, 22], the numerical evolutions would probably 
require different initial conditions to incorporate both the expanding and contracting parts of the wave-packets at 
the initial instance in "time" . This issue is solved in the isotropic case by including a scalar field that plays the role 
of a global internal clock. We thus feel that our results are not in conflict with [21, 22], and that all of our works 
requires further development of the quantum dynamics including the specification of a good internal clock, physical 
inner product, and Dirac observables. 

Additionally, the validity of the Kantowski-Sachs description of the interior could be questioned. Boundary con- 
ditions must be provided for the numerical evolution at the horizon which are prescribed using classical solutions. 
However, this could be modified when considering a true quantization of the full spacetimc. In particular, the im- 
proved quantization scheme described in section III B leads to modifications at the horizon which may require input 
from the full quantization. Additionally, by considering the homogeneous interior, we can not probe a realistic collapse 
scenario. 

Therefore, an interesting next step would be to apply the techniques described here to the inhomogeneous spherically 
symmetric quantization. The results described in [23] provide the initial steps at an amenable quantization, though 
the results presented there are only valid in the exterior. If the results can be extended to cover the whole spacetime, 
then an effective dynamics analysis could be extended. With this model, one could analyze a genuine collapse scenario. 
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APPENDIX A: ALTERNATIVE QUANTUM HAMILTONIAN 

Here we consider an alternate behavior for the S factors that are more heuristically motivated. In the language 
of [17] this corresponds to a lattice refinement model whereby the number of vertices is proportional to the transverse 
area. In addition this scheme has been considered in anisotropic Bianchi I models in [18, 19]. A stability analysis 
indicates that this scheme leads to an unstable difference equation [17] thus it may not represent a good quantization 
scheme. An advantage of this scheme is that the quantum difference equation is much simpler than in the previous 
one and more amenable to a proper analysis of the true quantum dynamics although the instability of the difference 
equation could spoil this. 
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In this scheme we have the 5 factors behaving as 



Again the effective Hamiltonian is given by 



Hg = - 



N 
2Ch 2 



5 b = 

So = 

sin bSh sin 5 r < 



A 

Pb 
A 
% 



fe+ ( 



sin 2 W 6 , 2 \ Pb 



+ 7 )—j= 



leading to equations of motion with lapse N = "f^/p^Sb/ (sin? 



c cos c5 c — 

2p c cos c5 c 
3 sin b5b 



4,6b 



1 



-p b cos ( 



sin c<5 c 

b cos 6(5& 

4 4 sin M h 

7 2 p b 6% cos b(5 b 

2 sin 2 bSh 



j 2 5b 7 2 &J 2 cos6<5b 



4 sin^ b5 b 



(Al) 
(A2) 

(A3) 



(A4) 
(A5) 
(A6) 

(A7) 



Again we can examine where the quantum corrections are appreciable by determining where b5b and c5 c are on 
the order of one. From the classical behavior, both 6<5b and c6 c blow up at the singularity and thus quantum effects 
are expected there. At the horizon, b5b vanishes and c5 c is small provided A is small. Hence with this quantization 
scheme, the effects near the horizon are minimal. The resulting numerical solutions are qualitatively very similar to 
those obtained for the constant 5 quantization scheme. As discussed above, the parameter controls the mass of 
the second black hole in the effective quantum theory which can be tuned to give a symmetric spacetime. We note 
that similarly to the other quantization schemes, the singularity is resolved dynamically. 




FIG. 7: These figures show the effective plots of pt = pb(p c )- The dashed lines represent the classical solution while the 
solid lines correspond to the solution of the effective quantum theory. The initial singularity is removed and the spacetime 
is connected to a second black hole. The three plots on the left show the solution for varying mass parameter, m = 1; 1.5; 2 
respectively. On the right p£ is varied in the three plots, p^ = 1; 2/3; 1/3 respectively. 
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